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Mapping out of equilibrium into equilibrium in one-dimensional 
transport models. 

Abstract. Systems with conserved currents driven by reservoirs at the boundaries offer an opportunity 
for a general analytic study that is unparalleled in more general out of equilibrium systems. The evolution 
of coarse-grained variables is governed by stochastic hydrodynamic equations in the limit of small noise. 
As such it is amenable to a treatment formally equal to the semiclassical limit of quantum mechanics, 
which reduces the problem of finding the full distribution functions to the solution of a set of Hamiltonian 
equations. It is in general not possible to solve such equations explicitly, but for an interesting set of problems 
(driven Symmetric Exclusion Process and Kipnis-Marchioro-Presutti model) it can be done by a sequence 
of remarkable changes of variables. We show that at the bottom of this 'miracle' is the surprising fact 
that these models can be taken through a non-local transformation into isolated systems satisfying detailed 
balance, with probability distribution given by the Gibbs-Boltzmann measure. This procedure can in fact 
also be used to obtain an elegant solution of the much simpler problem of non-interacting particles diffusing 
in a one-dimensional potential, again using a transformation that maps the driven problem into an undriven 
one. 



PACS numbers: 02.70.-c, 05.70.Ln 



1. Introduction 

Transport models are systems with conserved currents. In certain cases, they are such that their evolution 
leads to equilibrium when they arc isolated or in contact with a thermal bath. The probability distribution is 
then of the Gibbs-Boltzmann form. Coupling the boundaries to several external sources may induce currents 
across the bulk, driving the system out of equilibrium. In that case we do not have any general explicit 
formula for the distribution of probability of configurations, even in the stationary regime reached after long 
times. 

In order to make progress, one strategy has been to study systems that, due to their specific symmetries, 
admit a complete solution. Thus, in recent years a number of remarkable analytic results have been found for 
simple transport models (see [11] and references therein). Two important examples are the Simple Symmetric 
Exclusion Model (SSEP) — a one-dimensional system of particles, and the Kipnis-Marchioro-Presutti model 
(KMP) [27], a model of energy transport. For the former, Derrida, Lebowitz and Speer [10] (DLS) obtained 
an exact expression for the large deviation function of density using a matrix method that had been developed 
previously [9]. 

An alternative strategy is to restrict the calculation to the probability distributions of coarse-grained 
variables. If one considers conserved quantities, then the macroscopic fluctuations obey hydrodynamic 
equations with noise, the latter a manifestation of the microscopic chaos or stochasticity. Clearly, the 
more coarse-grained the description, the lower the level of noise, because fluctuations tend to average away. 
One is thus lead, in the macroscopic limit, to deterministic equations perturbed by stochastic terms whose 
variance is of the order of the inverse of coarse-graining box size N. As usual, one can recast the problem in 
terms of the evolution in time of the probability distribution. This is given by the Fokker-Planck equation, 
which is closely analogous to a Schrodinger equation in imaginary time, with the small parameter N^ 1 
playing the role of h. The 'semiclassical' treatment of these equations [16] follows the same lines as the 
derivation of classical from quantum mechanics, or geometric optics from wave dynamics. The logarithm of 
the transition probability obeys a Hamilton-Jacobi equation whose characteristics are trajectories satisfying 
Hamilton's equations. Macroscopic Fluctuation theory, developed by Bertini, De Sole, Gabrielli, Jona-Lasinio 
and Landim (BDGJL) [1, 3] is the resulting classical Hamiltonian field theory describing coarse-grained 
fluctuations. 

Up to this point the formalism is completely general. However, an analytic expression for the solutions 
of the classical equations is not possible for every model, so that even in the coarse-grained limit the problem 
has no closed solution. Remarkably for the hydrodynamic limit of the driven SSEP [1], BDGJL were 
able to integrate explicitly the corresponding Hamilton-Jacobi equations and recover the large-deviation 
function. They thus followed a path that is in principle logically independent of the one used to obtain the 
exact microscopic solution. Their derivation amounts to rewriting the problem in a carefully chosen set of 
variables. An analogous strategy subsequently allowed Bertini, Gabrielli and Lebowitz [3] to do the same 
for the hydrodynamic limit of KMP. 

One is left wondering what is the underlying reason for the existence of changes of variable that allow to 
completely solve the Hamiltonian equations, and how general their applicability is. In this paper we show that 
in the cases where this has been possible, there exists a non-local mapping taking the hydrodynamic equations 
of the model in contact with reservoirs into those of an isolated, equilibrium system. Large deviations and 
optimal trajectories are easily obtained in this representation using the detailed balance property, and can 
then be mapped back to the original setting in which detailed balance was broken by the boundary conditions. 
In the transformed, isolated model, spontaneous rare fluctuations are the time-reverse of relaxations to the 
average profile, but this symmetry is lost (as it should) in the mapping back to the original model. This 
accounts in this case for breaking of the Onsager-Machlup symmetry [31] between birth and death of a 
fluctuation, which has received considerable interest [29, 6] over the past few years. A short account of this 
work has appeared elsewhere [38] . 

The layout of this article is as follows. In section 2 we review the expression of particle exclusion problem 
in terms of spin operators. We construct a (coherent-state) path integral for the transition probability and 
derive from it the hydrodynamic limit. This route to the hydrodynamic limit is conceptually simple, and in 
addition has the advantage that the spin notation makes the symmetries and integrability properties of the 
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hydrodynamic limit explicit. Perhaps more surprisingly, the coherent-state representation yields directly a 
set of (Doi-Peliti) variables in terms of which the problem can be explicitly solved. 

In section 3 we take the more direct route starting directly from Fluctuating Hydrodynamics, and 
deriving from it the coarse-grained hydrodynamic equations. The reader may skip section 2 and start here, 
the only loss is that the symmetries of the hydrodynamic equations are not explicit and the boundary 
conditions less straightforward. 

In section 4 we review how the original stochastic dynamics for a single density field p{x) leads, in the 
low-noise limit, to a classical Hamiltonian field theory in a phase-space with two fields p(x) and p(x). The 
large deviation function - the logarithm of the probability of a configuration in the stationary regime - is 
given by the action of a 'classical' trajectory \p(x,t), p{x,t)] reaching the configuration at very long times. 
The formalism is a variant of the low-noise Frcidlin-Wcntzcll formalism, itself an implementation of the usual 
WKB semiclassical theory. 

For an undriven system with Detailed Balance (section 5) the dynamics has a symmetry between 
'downhill' relaxations and 'uphill' excursions. Using this symmetry one obtains the uphill trajectories as 
the time-reversed of the downhill ones, and this allows to compute explicitly the large-deviation function. In 
the case of driven transport models, the boundary terms violate detailed balance, and there is no obvious 
symmetry playing the role of time-reversal, no general method to find the uphill trajectories, and hence no 
explicit solution for the large-deviation function. 

In section 6, a paraphrase of the solution of BDGJL, we show that in the particular case of the SSEP, 
there is a very special set of variables that allows to solve completely the driven problem. The main result of 
this paper, in section 7, is to show that at the bottom of this possibility is the fact that there is a non-local 
mapping converting the driven chain into an isolated, undriven one. In section 8 we briefly show that one 
can apply the same arguments to solve the KMP model. 

To conclude, in section 9 we discuss how the method can be applied to the much simpler case of 
non-interacting particles diffusing in a generic potential, driven out-of-equilibrium by boundary terms. 

1.1. Mapping a driven into an undriven problem 

Let us consider here the simplest example of a mapping from out of equilibrium to equilibrium. This exercise 
will help to fix ideas, and to see how much one gets out of such an approach. 

Consider a diffusion process given by a Fokker-Planck equation driven by boundary conditions P(0) = p 
and P{L) = pl 

P- d ( T d + dVS ] P (1) 
dx \ dx dx ) 

Because of the boundary conditions, the current J = {—T-^ — 4^£) P is in general non-zero. 

Now let us introduce e@ v P = P\. Equation (1) maps to the backward Fokker-Planck equation: 

P - (t d dV \ d P (2) 
\ dx dx J dx 

Defining P> = £P lt we get: 

P' - d (t d dV ^\ P' (3) 

dx \ dx dx J 

This takes the form of a probability P' describing an evolution in a potential —V. The remarkable fact is 
that the original boundary condition on P expresses in the new variables just a normalization for P'\ 

'' L P'dx = e^pL - e^^po (4) 



Furthermore, it also implies that the current J' associated with P' vanishes at the ends A = (0,L), since 

P(A) = - ± ( T ± + P(A) = - ( T ± d f) P'\ A , = -./'( A i = o 
dx \ dx dx I \ dx dx 
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P' thus evolves with a Fokker-Planck equation with potential —V and no current at the boundaries; the integral 
of P' is conserved and the process satisfies detailed balance. The stationary measure is thus P' stat (x) oc e^ v ^ x > 

P' tat oc - P stat (x) = Poe-^- V ^ ] + c f e^^^dx' (6) 



Jo 

where the constant c is fixed by the right boundary condition P(L) = pi,. The overall distribution then 
reads 

_ p oe -/3[V(x)-V(0)] JL e fW( X ') dx , + pLe -0[V{ X )-V{L)) J* e PV( X ') dx , 

Pstatyx) — , (7) 

j^ e PV(x') dx f 

This result could have been obtained from the beginning by quadratures [25]. Actually, we have obtained 
much more, since we have mapped the evolution operator into the one of an equilibrium problem and we now 
understand the time-evolution of the primed variable as the relaxation of an isolated system: even without 
calculating anything, we have an intuition of the qualitative behavior. 

To see that the transformation is non-local, we may consider the expectation value of a function 0(x): 

(0)(t)= / dxO{x)P{x)= dxO{x)e- f3vi - x) P 1 {x) 
Jo Jo 

L 

dxO'(x)P'(x) - 0'{L)p L e fw(L) (8) 

All the time dependence is given by the expectation value of the new operator: 

0'{x) = - [ dy 0(y)e~^ VM (9) 
Jo 

which is a non-local function of the original one. 

In section 9 we shall show that actually one can use the same non-local transformation to map any single- 
particle diffusion model in a one-dimensional potential with sources at the ends into the same equilibrium 
diffusion problem with no sources. 

In the rest of the paper, we shall meet an analogous situation, for models with many interacting particles. 
The main difference is that we shall be making a transformation at the level of fields, not of their probability 
(P here is a probability of the O-dimensional 'field' x) . 
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2. Path-integral representations of exclusion processes 

In this section we introduce the exclusion processes. We write their evolution matrix in terms of (quantum) 
spin operators. This makes explicit the fact that the problem has more symmetries than mere particle 
conservation. Using standard spin coherent state techniques we present a novel derivation of the path 
integral and then take the hydrodynamic limit. We also obtain a natural set of variables F, F (related to 
the stereographic representation of the spins) which are often implicitly used in the literature. Finally, we 
write the hydrodynamic equations in terms of the average particle density p, and its conjugate variable p. 

The main result of this section is the construction of an hydrodynamic action which gives the logarithm of 
the transition probability between two smooth density profiles, together with the associated spatio-temporal 
boundary conditions. It is given by equations (46-49). The reader familiar with these representations of 
exclusion processes may skip this section, and find the more standard approach in the next one. 



2.1. Partial Exclusion Process 

Let us consider the symmetric partial exclusion process, a generalization of simple symmetric exclusion 
processes introduced in [35]. It consists of a one dimensional lattice gas, for which all sites can be occupied 
by at most 2j particles. The probability of a jump between a site and its neighbor is proportional to both 
the occupation number of the starting site and the proportion of vacancies of the target one]:: 

W(. . . , n k ,n k +i, n k + 1, n k+1 - 1, . . .) = — (2j - n k ) n k+1 (10) 

2j 

P 

W(. . . , n k , n k+ i, n fc - 1, n k+1 + 1, .,.) = — n k (2j - n k+ \) 

We fix a time-scale by choosing p — 1/2. The system can be put in contact at sites 1 and L with reservoirs 
of densities po and p\. This is usually done by introducing four rates a, 6 and 7, /3 which correspond to 
deposition and evaporation of particles at site 1 and L, respectively. We thus have the added rates of 
interchange with the reservoirs 

W{n u . . . -► m + 1, . . .) = a(2j - m) W(. . . ,n L -» . . . ,n L + 1) = 6 (2j - n L ) 
W(rii,... ->• m - 1,...) = 7m W(. ..,n L -> ...,n L - 1) = j3n L (11) 

Though the bulk-diffusion is symmetric, and the system thus satisfies a local detail balance relation, it 
can be driven out of equilibrium by the boundaries, if the densities imposed by the reservoirs are different: 

ol , 8 ,„„. 

Po = —— *Pi = j—, 12 

a + 7 + p 

These models are amongst the simplest interacting many particle systems driven out of equilibrium by the 
boundary sources. A schematic representation of the partial exclusion process is shown on figure 1. For 
j = 1/2, we recover the usual SSEP, which is known to be related to the 1/2 representation of the SU(2) 
group, whereas the partial exclusion processes correspond to the spin j representation [35]. For sake of 
completeness, we give the details of the relations with the spin operators in the next subsections and use the 
SU(2) coherent states to construct a path-integral representation afterwards. 



2.2. Master equation and spin representation 

The evolution of the probability P{n) of observing a configuration defined by the vector of occupation 
numbers n = (rii, . . . ,til) is given by the master equation 

= £ W ( n ' n)P(n') - W{n -> n')P(n) (13) 

Tl.' 7^71 

| The prefactor l/(2j) ensures that the rate at which a fixed number of particles jump to an empty site does not diverge with 
j- 
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Figure 1. Schematic representation of a partial exclusion process for j=3/2. There can be at most 3 
particles per site. Particles are injected at site 1 and L with rate a and <5 and can jump to the corresponding 
reservoirs with rate <5 and j3. 



where W(n' — > n) is the transition rate from configuration n' to configuration n. To keep the notation as 
compact as possible, we introduce = rife + 1 and = rife — 1. From (10) and (11), the master equation 
reads 

BPln) - ■ 



dt 



")n+ +1 P(. . . , n fc , n+ +1 , . . .) + n+(2j - n fc+1 )P(. . . , n+ n fc+1 , 



fc=i 



•■■) 

(14) 



- [(2j - n k )n k+ i + n k (2j - n k+ i)]P(. . . , n k ,n k+ i, . . .) 

+ a(2j - nr)P(nr, . . .) + jn+P(n+, . . .) + 5(2j - n£)P(. ..,n£) + /3n+P(. . . ,n+) 

- [7m + a(2j - m) + 6(2 j - n L ) + /3n L ]P(m,. . . ,n L ) 

The two first lines correspond to the dynamics in the bulk whereas the last ones stand for the interaction 
with the reservoirs. 

Let us now introduce a non-Hermitian representation of the SU(2) group to write the master equation 
in an operatorial form. A configuration of the system can be written as the tensor product of states of each 
sites i: \ip) = ®i\ipi), where each state \ipi) is given by a 2j + 1 components vector, such that an occupation 
number equal to n is represented by 



(15) 



where the 1 is on the (71 + l)th line, starting from the top. One then defines the (2j + 1) x (2j + 1) matrices: 
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Vo 
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• 3) 



(16) 



whose action on the state |n») are given by 
S+\ ni ) = (2j-m)\ni + l) 



S i \m) = rn\ni - 1) Sf\rii) = (n* - j)\n 2 ) 



Direct computations show that they satisfy the commutation relations 

[S^S^ =±S ± [S + ,S~] = 2S Z [S\S j ] =ie ilk S k Vi,j,k € {x,y,z} 



(17) 



(18) 



G 



They thus form a 2j + 1 dimensional representation of the SU(2) group, that is a representation of spin j, 
the usual magnetic number rm being related to the occupation number rij through m.j = rij — j. This is a 
non-unitary representation as «S + is not the adjoint of S~ §. 

The evolution operator of the partial exclusion process can now be written using these spin operators. 
To make this explicit, let us introduce the vector \tp) — J2 n P( n )\ n ) where the sum runs over all the possible 
configurations n — (n%, . . . ,tie,). Using the master equation, the time derivative of \ip) is given in terms of 
the spin operators by 





T~4 




H = H 1 +H B 






1 

= -Y\- 

4? ^ v 
■' fc=i 


&k- 


hl S+ + (j + S* k+1 )(j 


- s z k )] + [Sk+i s k + (i - s* k+1 )(j + si 


Hi 


= -a [Sf - 


U 


- Sf)} - 7 [ST - (j - 


+■ ^ 2 )] ; 


H l 


= -s [sl - 




- SI)] -f3[SZ- (j - 


+■ ^£)] • 



(19) 

-Hb corresponds to the dynamics in the bulk, whereas H\ and Hl result from the coupling with the reservoirs 
on the sites 1 and L. Introducing spin vectors Sk = {S^, Si, S%), Hb can be written in a more compact way 
as 

HB = --Y,(Sk-S k+1 -j 2 ), (20) 

3 i=l 

This is the usual connexion between exclusion processes and spin chains [12, 18, 35, 14]. 



2.3. Path integral representation in the hydrodynamic limit 

The hydrodynamic limit of lattice-gas models is usually achieved through the definition of a coarse-grained 
density field p, averaged over macroscopic boxes, the size of which is then sent to infinity [37]. Here, each 
site of the lattice already contains up to 2j particles. We can thus obtain a hydrodynamic limit by taking 
the limit of j, L going to infinity. In the macroscopic limit, all the spin representations are equivalent, j and 
L enter through the combination jL. This is just the manifestation of the fact [14] that the hydrodynamic 
variables represent the total spin in a box, whatever the spin of the elementary sites. 

2.3.1. Coherent states. To construct a path integral representation of H, we shall use, for each site k of 
the lattice, the following right and left spin coherent states [33] 



For an extended system of k sites, one introduces the tensor product 

l*)=®l**>. (22) 

k 

Due to the non-Hermiticity of the representation (16), \zk) is not the adjoint of (zk\. The construction 
of the path integral relies on the following representation of the identity 

dfi(z k ) \z k )(z k \ = 1 with dfi(z k ) = 2j ' + 1 d Zk _ , (23) 

7T (i + z k z k y 

§ Note that a similar path could be followed with the usual unitary representation of SU(2), see Appendix C. It would however 
lead to additional temporal boundary terms in the path integral, which makes the analysis of the action less straightforward. 
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2.3.2. Action functional in the large spin limit. The functional approach for spin operators has some 
subtleties [17] (see Solari [36], Kochetov [28], Vieira and Sacramento [39] for a derivation of the action and 
of the associated time boundary conditions) . For sake of clarity, most of the technical details are presented 
in Appendix A and we simply outline below the main steps in the path-integral construction. 

Keeping in mind that we ultimately want to describe a continuum theory, we replace number occupations 
by 'discrete' densities: 



Pk =2j 



(24) 



We would like to compute P(pf ,T; p 1 , 0), the probability of observing the system in state (p{, . . . , pi) 
at time T, starting from (p\, . . . , p % L ) at time 0, that is the propagator between two states (p* | and \p' 1 ) with 
fixed initial and final number of particles in each site. Using 2L representations of the identity (23), we write 



P(pf,T;p\0) 



(pV T *|p«) 



ndM4)dM(4)<PV><zV T V><zV 



(25) 



Because we are interested in the hydrodynamic limit (large jL), we may first take a large j limit and then 
send L to oo. In this 'large spin' limit, (25) reads (see details in Appendix A) 



P(pf,T;p\0)= / JJdM(4)d M (4) / VzVz U S [i 



ZkZk 

1 + z k z k 



1 / 



logZfc - l0g(l + ZfeZfe) 



e=i,/ 



2 j / dt 



p% exp[-5], 



(26) 



El 



ZkZk 



ZkZk 



H(z,z) 



(27) 



The role of the Hamiltonian H(z, z) is played by the quantity —-^{z\H\z) which is computed using [17]: 



{zk\St \zk) = 2j- 



Zk 



(z k \S k \z k ) = 2j- 



z k 



k ' J l + z k z k 

Explicitly, this yields 

H(z, z) = Hb{z 7 z) + Ho(z, z) + Hl{z 7 z) 

z k Z k+ l {zk+i — z k ) 2 



z k z k 



{z k \S z k \z k ) = j 



z k z k - 1 
z k z k + 1 



L-l 

H B {z,z) = — 



H {zi,zi) = a 



k=l 



+- z k z k 1 + z k+ iz k+ i 2 

1 Zl(l - Zl) 

+ 7 ; , : _ ; n L {zL,ZL) 



Z k +1 



2fc 



l + Z k+1 Z k+1 l + Z k Z k 



ZL 



1 



1 + ZlZl 1 + ZlZl ' ' " l + ZLZL l + ZLZL 

Initial and final conditions on the fields z, z are imposed by the delta functions in (26) and read 



1 



Pk 



-f f 

44 



(28) 

(29) 

Zk+l — z k 
2 

(30) 

(31) 



2.3.3. Density field. To get more insight on the physics of this field theory, we introduce a new 
parametrization 



so that 



Zk 



p k 



1 - Pk 



Zk 



ZkZk 



z k z k 2j 



(32) 



(33) 



plays the role of a density. The transformation (32) is such that p is canonically conjugated to p. The time 
boundary conditions (31) on the field can be written, as expected, as 

Pk(0)=pi, Pk(T) = p{, (34) 
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whereas /3fc(0) and p k (T) are unconstrained (for details, see Appendix A). This highlights the correspondence 
between the field p and the actual density of the system, as do (33). From (28) and (32), one sees the 
correspondence between spins and densities 

(z k \S+\zk)=2j(l-p k )e^, (z k \S^\z k ) = 2jp k e-P k , (z k \S z k \z k ) = j(2p k - 1)(35) 

This yields for the Hamiltonian 

1 

H (P> P) = o E U 1 - ^ W [e h ~ pk+1 - 1] + p fc (l - Pfe+i) [e^ 1 -^ - 1] } (36) 



fc=i 



+ {a(l - pi)(e^ - 1) + 7pi(e~* - 1) + <5(1 - Pi )(e* £ - 1) + ^(e"^ - 1)} 
Making the change of variables in the action gives 

P(n f ,T;m,0) = JvpVpexp{~S[p,p}} 



S[p,p] 



2j 



df 



E PkPk - W(p, p) 



where one integrates over fields p satisfying the temporal boundary conditions (34). 



(37) 
(38) 



2.3.4- Connection with Doi-Peliti variables In the previous subsections we took advantage of the SU(2) 
coherent states to construct the path-integral representation of the evolution operator. A more standard 
approach, based on Doi-Peliti formalism, could have been followed. Starting from the density fields p, p, the 
usual bosonic coherent states can be obtained through a Cole-Hopf transformation [4] : 

P = U*; P = log<f (39) 
which in the z, z variables reads 

cj>* = z; 4>^t^-. (40) 
1 + zz 

From (29) one sees that the bulk Hamiltonian Hb reads in these variables 

n B (cp, 4>*) - — ^Mk+ii^U! - & + (4> k +i - - <t>l) (4i) 

k 

To understand the meaning of the <j), <f>* variables, we recall the action of the spin operators in the 
coherent state representation: 



(*k\s£\4>) = 



2 d 

2j ~ Zk -~ Zk W k 



dz k 



z k 



d 

dz k 



3 



(42) 



(44) 



The matrix element used to construct the path integral (28) are given in terms of (/>, (jf by 

(z k \S+\z k ) = 2j<f>l-{cj>i) 2 {2j4> k ) ] (z k \S k \z k ) = 2j<j> k ; (z k \S z k \z k ) = <f>* k (2j<l> k )-j , (43) 
and one sees by comparing (42) and (43) that 

4> and <j>* satisfy the usual bosonic commutation relation. They correspond to the usual Doi-Peliti [13, 32] 
representation of bosons. We thus see that one can construct a path integral in terms of the variables <f> and 
<f>* either by a Cole-Hopf transformation, or by directly expressing the spin operators in the representation 
(42) and then constructing the path integral treating z and as conjugate bosons ||. In what follows, 

however, it will be useful not to lose sight of the SU(2) symmetry. 

|| One must however proceed with care as states with more than 2j particles are not physical and are difficult to handle using 
bosonic coherent state. See for instance [40]. 
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L — l i 

. 7 E- dx ( 45 ) 



2.3.5. Hydrodynamic limit So far, the action is still a discrete sum over the whole lattice. We shall now 
take the full hydrodynamic limit to describe the evolution of smooth profiles on diffusive timescales. We 
thus introduce a space parametrization Xk = j; and rescale the time t — > l?t. At the macroscopic level, the 
density profiles are smooth functions and discrete gradients can be replaced by continuous ones: 

Vp(x fe ) . . Vp{x k ) ' 1 1 

Pk+i - pk —> — 7 — , pk+i - Pk — * — ^ — , L 

k=l " u 

For a more rigorous approach, see [37, 2]. The first order in a j; expansion of the action then reads 

P(p f ,tf,p i ,0) = JvpVp C - 2: > LS[p P ] (46) 

S[p,p] = J ' j dxdt{pd t p~Hlp,p}}; n[p,p} = ^a(S7p) 2 -^ P S7p (47) 
where a = p(l — p). 

In Appendix B we show that fields p and p are constrained, in the hydrodynamic limit, to satisfy the 
spatio-temporal boundary conditions 

a 6 

Vt, p(0,<)=p -^— , p(M =Pi = JXfl' P(°)=P(L =0 (48) 
a + 7 o + p 

p(x,0) = Pi (x), p(x,T)=p f (x) (49) 

If a,f3,j,S ~ 0(1) at the microscopic level, these conditions are strict, in the sense that the fields do not 
fluctuate in the borders ^f, whereas a, /3,7, S ~ ©(L^ 1 ) would also allow fluctuations of p and p at the 
boundaries. 

Due to the correspondence (20), (35) with the spin operators, the whole process described here simply 
amounts to taking the classical limit of a Heisenberg spin chain with some particular boundary conditions. 
One can indeed check that the Hamiltonian (47) corresponds to 

Ub = -\ [ dxVS- (50) 



where the classical spins [45] are defined as S x ' v,z = linx 7 -_ > . 0O (2j) 1 (z\S x,y,z \z) . The spatial boundary 
conditions in terms of classical spins are then given by 

S z {Q)=Po-\ S z {\) = Pl - 1 - (51) 

S x (0) = \ S*(l) = ~ 

sm = S y {1) = ±Z*EL 

S+(0) = l-p o S+(l) = l- Pl 

S-(0)=p o S-(l)= Pl 

% The probability to observe a smooth profile such that p(0) ^ po scales as e~ jL . 
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3. Fluctuating hydrodynamics 

Let us briefly review in this section the construction of the action starting from fluctuating hydrodynamics 
and using the Martin-Siggia-Rose [30] , DeDominicis- Janssen [7, 22] formalism. 

In terms of the instantaneous current J(x) at site x, defined by the continuity equation, the evolution 
of the density is given by the stochastic equation 

P=-VJ; J = -\v P -^ m p(0) = p ; p{l)=Pi (52) 

where 7/ is a white noise of variance 1/ (2j L) and a = p(l — p). This is the usual formula for the fluctuating 
hydrodynamics of the exclusion process [37] . We write this as a sum over paths and noise realizations with 
a delta function imposing the equations of motion and a Gaussian weight for the noise: 

P{p f ,t f ;p\ 0) = J VrfDp 6[p + V J] e~ 2j L & ft dxdt (53) 

Exponentiating the functional delta function with the aid of a function p(x,t): 

where p is integrated along the imaginary axis. Integrating by parts, we obtain: 

P(pf,tf,p\0) = J Vr]VpVpe- 2jL ^ ft*** {pp+Vp{^p+V^v))+W{^)} (55) 

We can now integrate away the noise: 

P(pf,t f ;p i ,a) = JvpVpe- 2 ^ S^dxdtipp+ivpvp-iivp)^} (56) 

which reads 

P{p s ,t s ;p\0) = JvpVp e - 2jL V° f ti dxdt (pp- n » (57) 

to obtain 

H = \ [a(V/5) 2 - VpVp] (58) 

which is equivalent to (4G). The paths are constrained to be Pi{x) and Pf{x) at initial and final times, 
respectively. The values of p are unconstrained, which is in agreement with the fact that this is a Hamiltonian 
problem with two sets of boundary conditions. The construction above can thus be seen as a formal Hubbard- 
Stratonovich transformation to introduce the p field. 

Let us finally note that from equation (52) and (56), one sees that 

P(pf,t f ;p\0) = y*2 3/5 2 3j5 e- 2 ^/o / /o 1 dxdt(JVp+iVpVp-i(Vp)V} (5Q) 
Formally integrating over V/5 gives back the usual fluctuating hydrodynamics [1, 24, 34] 

P(pf,t f ;p\0)= J Vp exp [ - 2jl£ j dxdt (J + ^ /2)2 ] (60) 

In all this we have been very sloppy about the spatial conditions p should satisfy: see Appendix B.2 for 
details. 
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4. Large deviations in the coarse-grained limit 



In this section we review the steps leading from the fluctuating theory to a non-fluctuating Hamiltonian 
dynamics, valid in the low-noise limit — itself arising in the large coarse- graining limit. 

4-.1. Classical solutions 

In order to calculate a probability we have to evaluate 

P(pf, T)= J dp'Pipf^T; p\ 0)P(p\ 0) (61) 

which is true for any time T. Using the path integral expressions of the previous sections (cfr (26)), one 
then gets the sum of the exponential of the action 

P(p f ,T) = J dp 1 J V[p,p]e^ LS ^P{p\0) (62) 

over trajectories with initial and final profiles p 1 and pf . To leading order in jL, we have that P{p? , T: p l , 0) 
is dominated by the trajectories extremalizing the action (46) 

S = J dxdt [pp - H] (63) 

i.e. satisfying Hamilton's equations: 

6 f 

P(x, t) = j^— / dx'dt"H[p(x', t'),p(x', t')} (64) 
op[x,t) J 

P&t) = - Tl ^ / dx'dt'H[p[x',^),p{x',if)} (65) 



5p(x,t) 

These are completely determined (at least up to a discrete set of trajectories) by the initial and final values 
of p{x). What we have outlined is the exact analogue of the way classical trajectories dominate the path- 
integral in the semi-classical limit h — > in quantum mechanics. Similar approaches have been used many 
times, as for instance to analyze the noisy Burgers equation [15] or in reaction diffusion systems [26]. 
In the case of the SSEP, Hamilton's equations read: 

P (x) = ^Ap-W[aWp] (66) 

'p{x) = {p-\){Vpf-\llp (67) 

(66) is a conservation equation and thus defines a current through p = — V J p where 

J p \p,p] = ~Vp + aVp (68) 

To determine the probability of a transition between a profile p z and another profile p-" in a time T, one thus 
has to find the trajectory (p(x, t), p(x, t)) that solves (66) and satisfies the appropriate boundary conditions. 
The action of this trajectory then yields the logarithm of the probability of the transition. 

We are thus led to solving a classical Hamiltonian field problem, where initial and final positions are 
fixed and momenta unknown. This is in general very difficult - even numerically, where one has to solve a 
'shooting' problem to reach the desired final configuration at the right time [47]. 

For completeness, let us write the equations of motion in terms of classical spins satisfying the Poisson 
bracket algebra [46]: 

{S z (x),S ± (y)}=±S ± (x)6(x-y) ; {S+(x), S~(y)} = 2S Z 5(x - y) (69) 

which corresponds to the quantum commutations (18) and yields the equation of motion in terms of spin 
variables through 

S = {S,H} = iASAS (70) 
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4-2. Downhill trajectories 



For classical equations deriving from a stochastic problem, there always exists a class of trajectories which 
are easy to find: those that arc overwhelmingly the most likely in the low noise limit. Here and in what 
follows, we shall call these 'downhill' trajectories. If one remembers that the dynamical action corresponds 
to the stochastic equation (52), one obtains these solutions by directly putting rj = 0. They correspond to 
the following solutions of the system (66): 

p{x) = \^P ftx,t) = (71) 

The solution corresponds to diffusive relaxation towards the linear stationary profile p(x) = (1 — x)po + xp\. 
The corresponding action is zero, in agreement with the fact that the corresponding probability is 1, which 
simply means that an initial configuration p(x) almost surely relaxes towards the stationary state. The 
stationary profile p = 0, p = p is a hyperbolic fixed point of the dynamics in the full phase-space, since it is 
missed as soon as p ^= in the initial condition. 



4- 3. Large deviation function from extremal trajectories 

As shown by equation (61), the probability P{p*) ~ e~ N:F ^ p ' to observe a profile p* is the average probability 
of going from an initial profile p i to the profile p* . The logarithm of this transition probability is given in 
the large N limit by the 'classical' action of a trajectory starting in p 1 and arriving in a time T in the 
configuration pf = p* , and satisfying (66). 

How can a trajectory just reach a generic configuration p* at a very large time T — > oo? The only 
possibility is that it takes a hyperbolic trajectory that falls in a finite time in the vicinity of the stationary 
profile, stays there almost all the time, and then goes to the profile p* in a finite time. This is illustrated 
in figure 2: trajectories that matter at long times are thus near-misses of the stationary points. The first 
part of such a trajectory (essentially the diffusive relaxation towards the stationary profile) has almost zero 
action. Thus P(p z , — > p* , T) — P(p, p* , T') where T" is a time which differs from T by a finite contribution, 
not relevant in the T — > oo limit: the transition probability P(p l , — > p* , T) is independent of p 1 in the long 
time limit and equation (61) becomes 

P(p*,T^oo) = J dp i P(p i ,0)P(p,p* 1 T' - oo) = P(p,p*,T' oo) (72) 

To determine the probability to observe a given profile p* , one thus has to find the extremal trajectory 
which starts at t = — oo in the stationary profile and arrives at t — in the desired profile p* . Its action 
then yields the large deviation function J-[p*] 

P(p*) = PGM = -oo,pV = 0) (73) 

We have assumed here that there are no metastable states. In cases in which many such states exist, one has 
to consider trajectories falling into each one of them, and also making jumps between them before reaching 
the final point. 
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Figure 2. 'Downhill' (noiseless) trajectories have pit) = 0. Trajectories reaching an arbitrary point at 
long times are "near-misses" from the stationary point, and can be decomposed in the long-time limit in a 
'downhill', followed by an 'uphill' trajectory flowing into, and out of the stationary point, respectively. 



5. Detailed balance relation, relaxation-excursion symmetry 

Finding the trajectory that reaches a given density profile starting from the vicinity of the stationary one is 
in general a difficult task. There is however a class of systems for which this calculation greatly simplifies: 
those for which there is a detailed balance symmetry - playing the role of a time-reversal - that relates the 
paths followed by the system in a rare (noise-induced) excursion with the relaxation back to equilibrium. 
This Onsager-Machlup symmetry allows one to compute the rare excursions (in general difficult) from the 
sole knowledge of the relaxations (easy, as explained in the previous section). This is for instance the case 
for a SSEP in contact with reservoirs of equal densities po — pi, as we shall see below. 

At the level of operators, detailed balance simply says that the evolution operator and its adjoint are 
related by a similarity transformation. Here, we arc only interested in symmetries at the level of the action, 
which can be read as a canonical transformation followed by time reversal, leaving the action invariant up 
to boundary terms (See Appendix D). 

To make this explicit in our case, we write the Hamiltonian density as 

H[p, p]=\vp <rV [p - log j!L-) = \ Vp aV [p - (74) 
where V p is the equilibrium entropy: 

V p = j &x[p logp + (1 - p) log(l - p)} (75) 
The transformation we are looking for is given in two steps [23]: i) the canonical transformation 

p + ~r- ! P~^ P ( 76 ) 
dp 
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followed by: ii) a time reversal 

(p,p,t)^(-p,p,T-t) (77) 
In terms of the current (68), this reads: 

[p(x,t),J(x,t)} -» [p(x,-t),-J(x,-t)} (78) 
the meaning of which is transparent. The new, 'time-reversed' variables are: 

p TR (x,t) = p{x,-t) ; p TR (x,t) = -p(x,-t) + — -(a;,-t) (79) 

dp 

It is easy to check that (76,77) map the action into 

S[p, p] -> [lp]o + J dtdx {p TR p TR - H[p TR , p TR }} (80) 

The space and time boundary conditions are transformed from (48) and (49) to 

p TR (a;,0) = p*(x) p TR (x,T) = p(x) (81) 

p TR (0,t)=p p TIi (l,t) = px Ptr(0, t) = log - — — Atr(M) = log : Pl (82) 

1-/90 1 - pi 

The problem is thus recast into finding a trajectory starting in p* and relaxing towards the stationary profile. 

When the system is at equilibrium, in contact with two reservoirs imposing the same density on the two 

boundaries (po = pi), the zero-noise diffusive trajectory 

p TR (x, t) = C st = log p TR = i Ap TR (83) 

1 - Po 2 

is a legitimate solution of the classical equations which satisfies the new boundary conditions in space and 
time. The action of such a trajectory is 

dtdx {ptrPtr -H} = log P ° / dx(p - p*) (84) 



1 - Po 

Together with the boundary terms of (80), one gets for the large deviation function, in terms of the original 



varia 



iables 



S[p*} = / dx 



(1 -P) log^-^ + plog-^ 

1 - Po Po 



(85) 



which is the usual equilibrium entropy [2]. The computation of the extremal trajectory going from p to p* 
has thus been made possible by the time-reversal connection between excursion and relaxation induced by 
the detailed balance relation. Thanks to the mapping (76) and (77), one just has to find a trajectory going 
from p* to p - a relaxation - and from it obtain an excursion. 



5.1. System driven out of equilibrium - Violation of detailed balance and loss of Onsager-Machlup symmetry 

Let us now consider why this 'trick' does not work when the system is driven out-of-equilibrium by the 
boundaries. One can still make the transformation (76) and (77), as the bulk dynamics satisfies detailed 
balance. However, if p\ ^ po, (83) is not acceptable as p = (f does not satisfy the spatial boundary 
conditions. We conclude that the most probable excursion from p to p* is consequently not the time reversed 
of a diffusive relaxation: the Onsager-Machlup symmetry is broken by the boundaries. 
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6. Driven exclusion process: remarkable changes of variables 



As we have seen, in out of equilibrium systems detailed balance is violated - for example by the boundary 
conditions - and there is no simple relation between excursions and relaxations. One thus cannot, in general, 
compute easily the rare excursions away from the stationary state, and from it the large deviation function. 
Even if we have been able to reduce the computation of large-deviation functions to the solution of a problem 
of classical dynamics, we are not able to solve for its trajectories in a closed way. It may then come as a 
surprise that BDGJL [1, 2] were able to uncover what amounts to a scries of changes of variables, which 
end up by mapping the driven problem into one where there is a simple time-reversal symmetry between 
excursions and relaxations. This allowed them to compute the excursions in terms of the relaxations in the 
new variables, and then work their way back to the original variables in which the symmetry does not hold. 

In this section we shall paraphrase their derivation, to emphasize how surprising it is. In the next 
section we shall argue that at the bottom of this is the (also very surprising) fact that the SSEP driven out 
of equilibrium by the boundaries can be mapped back through a change of variables, in the hydrodynamic 
limit, to an equilibrium SSEP. 

Let us first note that the choice of axes we have used to write the hydrodynamic limit Hamiltonian (50) 
in terms of spin variables TLb — — h J dxVS* • V5 is arbitrary. Only the boundary conditions (51) break the 
rotation invariance. 

Let us make a transformation: 



S x — S z 



S' 



Sx 



s'„ 



-s. 



V 1 



(86) 



a rotation of angle 7r/2 around the y axis followed by a reflexion with respect to the x'-z' plane. In terms of 
the coherent-state coordinates - the stereographic representation of classical spins — this transformation is 
given by the simple homography 



z - 1 



z-1 



z + 1 z + 1 

In this set of variables, the action reads after a lengthy but straightforward computation 



(87) 



S= dx 



P lo gT ^ + (l-p)log^ 



S[z', 



2 2 

As this change of variable corresponds to a symmetry of the Hamiltonian, the action S[z' , z'\ is the same as 
the original one where z and z have been replaced by z' and z ' . It is thus obtained by taking the continuum 
limit of (26): 

S[z',z'} = j dxdt{J^ + \^zr(-^^\\vz^(-^—S\ (89) 



1 + z'z' 



Of course, we could have arrived at the same point by defining at the outset the coherent states in terms of 
the rotated operators. 

As previously, we can make a Cole-Hopf transformation to introduce Doi-Peliti like variables: 
1 + z z' 

Surprisingly, as we shall show below, this change of variables greatly simplifies the problem. To stay as close 
as possible to the solution introduced by BDGJL, we rather introduce a slightly different set of variables 

l + cb*' 



F = 



F = 2(f) 1 



(91) 



with which the same simplification occurs. In the F, F variables, the action is given by 



S = J dx 

S F [F,F] = I dtdx 



plog| + (l-p)log^-jl 



Sf[F,F] 



FF + ^F 2 WF 2 



-VFVF 



(92) 
(93) 
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The overall change of variables from p, p to F, F reads 

P 



F = 



p+(l-p)eP' 

In terms of these, the spatial boundary conditions read 

F{0)= Po F(l)= Pl F(0)=0 
while the equations of motion become 



F = (l-p){eP-l)- P (e-P-l) 
F(1)=0 



F = iAF-F(VF) 2 ; 



F = --AF- V 
2 



F 2 VF 



(94) 
(95) 
(96) 



In particular, by analogy with the 'downhill' zero-noise solutions of the previous sections, we can try 
F(t) = 0, which corresponds to 



F = -AF 
2 

In the original variables, this solution reads 



p = 0; 



(97) 



(98) 



which is nothing but the diffusive relaxation to the linearly stationary profile p. Diffusion in F, F variables 
thus corresponds to relaxation in p, p. 

6.1. A second detailed-balance like symmetry 

Remarkably, the action (92) has a detailed balance- like symmetry, which is unrelated to the original physical 
one. To see this, integrate the last term of Sf by parts [43], so that the action becomes: 



S F [F,F] = I dtdx { FF 
which can also be written 



1 



FVF 



S F [F,F] = / dtdx 



FF 



where we have introduced the potential 



1 



-FVF 



F 



F — 



AF 
VF 2 

SV F 

If 



dx logVF 



The composition of the transformations 



F 



p+ 5V F _ p 

f + Jf- f 



AF 



(F,t)^(-F,T-t) 



VF 2 ' 

is thus a symmetry of the action. Note the analogy with (76,77). The new variables read: 



F m (x,t) = F(x,-t) ; F TK (x,t) = -F(x,-t) + 



AF 

VF 2 



(x,-t) 



(99) 
(100) 
(101) 
(102) 
(103) 



Once again, rather than looking for excursions in the variables F, F, the problem is reduced to searching 
relaxations in the variable F TR , F TR . 

Because the system is driven out of equilibrium, it would be natural to expect, as in section 5, that 
this symmetry is broken. Remarkably, we shall show below that this symmetry is not violated by the spatial 
boundary conditions, in spite of the system being driven. Transformation (102) indeed maps the spatial 
boundary condition (95) into 



-F TR (0) = p ; F TR (l)=pi; 



F TR (0) 



AFr 



(VF TR ) 2 

Contrary to what happened in section (5), the 'zero noise solution' F 
satisfies (104). From (95), one indeed sees that 

F\ x =o,i = F x =o,i = 



AF TR 

= F ™ {1) -Wf^ 



= 

x = l 



(104) 



AF TR /2; F TR = this time 
(105) 
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Together with (96), this shows that any extremal trajectory satisfies 

AFU=o,i = 2[F + F(VF) 2 ]| X=(M = 
Furthermore, F TR = reads in the initial F, F variables 
AF 



(106) 

F -ww=° (107) 

which is indeed compatible with (105) and (106). Zero -noise diffusive relaxations in the variables F TR , F TR 
thus correspond to the excursion in the initial variables. The action of such a solution is 

T 

dadogVF (108) 



S F = [VfYo -- 

and the corresponding large deviation function is given, in the original variables, by 

n T 



S| "*' = / 



d.r 



p log | + (1 - p) log ^— ^ + log VF 



o 



F is determined from p by solving the equation F TR — which reads 

AF 



p = F + F{l-F) 



VF 2 



Injecting the temporal boundary conditions in (109), one gets 
S[p*] — dx 



P*log^ + (l-p*)log^-^ 



log- 



VF 



(109) 



(110) 



(111) 



Pi - Po_ 

This is indeed the solution of the problem, initially found by Derrida, Lebowitz and Speer [10] and later 
recovered by BDGJL [1]. Let us stress however that the existence of a potential functional Vf is a mystery, 
as nothing guarantees that such a function exists out of equilibrium. Yet another mystery is the fact that 
the symmetry related to Vf is unbroken by the reservoirs. We shall show below that these surprises are 
deeply related to the fact that this model can be brought back to equilibrium. 
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7. Non-local mapping to equilibrium 



In the previous section we have shown that if one rotates the axes, and expresses everything in terms of 
variables F, F associated with the spin operators in the coherent state representation (or, in the classical 
limit, with the stereographic projection of the spins), a new, surprising detailed-balance symmetry becomes 
explicit. It is unrelated to the original one and is not broken by the source terms at the boundaries. In 
this section we show that this sequence of miracles can be condensed into only one: at the hydrodynamic 
level the chain driven out of equilibrium can be mapped onto a free equilibrium chain with no sources at the 
boundaries. The transformation that does this is, however, non-local in space, and is the analogue of the 
one we discussed in the introduction for non-interacting particles in a potential. 
Starting from the action Sf (99), one introduces the non-local variables 

F' = \7F: \7F' = F A ' 



(VF) 2 

where the second equation can also be written: 



F 



F' 



1 

F' 



VF' 



VF' 



(112) 



(113) 



This change of variables takes the action into itself, apart from temporal boundary terms: 

T 

S= Idx 



plog | + (1 - p) log i— | + log F' - F'F' 



d.rdl <J F'F' + ^F' 2 (VF') 2 



-VF'X7F' 
2 



(114) 



Equation (112) thus describes a non-local symmetry of the Hamiltonian. This suggests that we continue the 
succession of changes of variables done up to here 

(p,p)^(F,F)^(F',F') (115) 

by a further transformation (F' , F') — » (p' , p'), where p' and p' are related to F' , F' in the same way as are 
the unprimed counterparts of (94) 



p > = F' + F'(1-F')F'; 
One thus ends up with 



p' = log(l+ ^ F 'A (116) 



1 - F'F' 



S= dx 



-S' 



plog 



(1 - p) log 



1 



F' 



1 - F 



1 - F' 



1 - F' 



logVf-^— ^-/Z log -(!-//) log ^—^ (117) 



(118) 

where S' = J dtdx {p' p' — H 1 } and Ti' is formally equivalent to the initial Hamiltonian density 7i but for the 
primed variables: 

1 _".2 1, 

2 

The overall change of variables, which reads 
1 



S' = / dtda; 



p'p' 



-a p ,Vp' +-Vp'Vp' 



(119) 



V 



1 - eP' 



(l-p)(ef-l)-p(e-"-l) 



(l-p')(e? - l)-p'{e-P - I) 



(120) 
(121) 



_p+{l-p)eP_ 

thus maps the action of the hydrodynamic limit of the SSEP into another SSEP. 

What we shall now prove is that this new process corresponds to an isolated chain, and consequently 
possesses a detailed balance relation, induced by its equilibrium entropy, which is not violated by the 
boundaries. Last, we shall show that this detailed balance relation, mapped back to the original p, p variables 
is the non-local symmetry (102). 
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7.1. Boundary conditions - Currents 



From (112) one sees that the spatial boundary conditions (95) read in the new variables 

= 0; / F' = Pl -p (122) 



F' 



In the language of spin variables, the hydrodynamic Hamiltonian (50) presents three conserved quantities in 
the bulk - the components of the spins - which can be written 

<2 1 =2p'-l = 2^ ; Q 2 = F' = 2(S' z +iS' y ) ; Q 3 = F'(l - 2F') = 2S' X - 1 (123) 

Their continuity equations read 

Qi = — VJi (124) 



where the currents Ji are given by 



J Ql = -Vp 1 + 2a p ,Vp' ■ J Q2 = l -VF' + F 2 VF' 

VP" - -2 

Jq 3 = (1 - IF') — + [F'+F' (1 - 2F')]VF' (125) 

Let us show that these currents vanish at the boundaries for all extremal trajectories. Such trajectories 
satisfy the equation of evolution (96) and the boundary condition (95). This implies that AF vanishes at 
the boundary (l.h.s. o f (96) together with F = 0). From the definition (112) one then gets that VF' also 
vanishes, which implies, together with the l.h.s. of (122), that VF' also vanishes. Last, from the mapping 
F',F' — > p' , p' one sees that if both VF' and VF' vanish, so do Vp' and Vp'. All in all, one gets that at the 
boundaries 

VF' = VF' = Vp' = Vp' = (126) 

We then see from (125) that the three currents Jq { vanish at the boundary: the transformed model in 
the primed variables is an isolated chain. This condition alone, supplemented with the r.h.s. of (122), 
encompasses all the original boundary conditions. 



7.2. Profiles and trajectories 

When the dust sets, we see that all that has been used is the fact that the original variables (p, p) or, more 
physically (p, J), have been mapped by a non-local transformation into new densities and currents 

{p(x,t),J(x,t)) -> {p'(x,t),J'(x,t)) (127) 

The original detailed balance symmetry (78) that maps a trajectory into its time-reversed: 

{p(x,t),J(x,t)) -> (p{x,-t),-J(x,-t)) (128) 

is broken by the source terms driving the system out of equilibrium. Miraculously, the transformed chain is 
isolated, so that time-reversed trajectories are related through: 

(p'(x,t),J'(x,t)) - (p'(x,-t),-J'(x,-t)) (129) 

Coming back to the original variables via (127) mixes the density p' (symmetric in time) with the current 
J' (antisymmetric in time), thus making the pair of transformed trajectories expressed in p and J neither 
symmetric nor antisymmetric. 

Let us now describe 'uphill' and 'downhill' trajectories. The stationary profile p maps to a flat profile 
-/ _ Qst ^ rp-^g p rec j se va i ue f (j st i s arbitrary, due to dilation- invariance of the model, contrary to F" which 
is constrained by the r.h.s of (122). 
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• Relaxations: diffusive trajectories of the initial model satisfy 

P=\^P, P = (130) 

From (120), one sees that p — implies V/5' = 0. As the primed variables also satisfy the equations of 
motion (66), the resulting trajectories evolve with 

p> = Up' (131) 

Relaxations thus map into relaxations. 

• Excursions: The instanton equations (107) imply VF' = (cfr (112)). Using the relation of F, F (94), 
as applied to the primed variables, one gets 

V/3' - ^ = (132) 
o- p , 

Injected back in the equations of motion (66), this shows that densities evolve with 

p' = --Ap' (133) 

Excursions of the initial model, once mapped back to equilibrium through (120), are given by time- 
reversal of the isolated chain's relaxations. The action S'[p',p') of such an uphill trajectory is 
J Q dx[p'logp' + (1 — //)log(l — p')]q ■ As F' is constant along the instanton and J Q p' is a constant of 
motion, the overall action (117) reduces to the large deviation function (111), as it should. 

7.3. Detailed balance symmetry 

Let us now show in terms of spin variables how the detailed balance relation of the isolated chain accounts 
for the miraculous transformation (102) which allowed BDGJL to compute the large deviation function. In 
the spin variables, the original detailed balance symmetry (76) and (77) amounts to a reflexion of all the 
spins with respect to the x — z plane: 

(S x ,S y ,S z )^(S x ,-S y ,S z ); T^T-t (134) 

Because the bulk Hamiltonian is also invariant with respect to any simultaneous rotation of all the spins, 
any composition of (134) with a rotation gives another 'detailed-balance like' symmetry. These symmetries 
are all broken by the boundary conditions of the original model. Once mapped to the isolated chain, the 
boundary conditions (122) reduce to fixing the value of an integral of motion: 

2 [ L (S' z +iS' v )= [ L F' = pi-p (135) 
Jo Jo 

Among all the 'detailed-balance like' symmetries of the isolated chain, only one preserves (135). We are thus 
left with the transformation 

{S x , S y , S z ) — > (—S x , S y , S z ); t — ► T — t (136) 
From the expression of Q3 in (123), one sees that 5 1 ™' = — S' x can be written as 

1 + F; R (1 - 20 = -1 - F"(l - 2F') (137) 
Using F' = VF and looking for Ftr = F, this reads 

F' + f; r = 1 + ^ (138) 
Differentiating once and using VF' = F — , one gets 

^* +#= (Hf (139) 

which is nothing but the non-local mapping (103) between 'downhill' diffusive solutions and the instantons. 
One thus sees that the miraculous solution of the initial model is simply induced by the detailed balance-like 
relation (136) of the isolated chain, which does not violate the boundary conditions (122). 
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8. KMP 



The Kipnis-Marchioro-Presutti model (KMP) was introduced in [27] as a one dimensional model of energy 
transport satisfying Fourier Law. Bertini, Gabrielli and Lebowitz recently computed the large deviation 
function of the energy profile using an approach similar to the one used previously for the SSEP [3]. We 
shall show below that once again a mapping back to equilibrium explains this success. 

The functional expression for the fluctuating hydrodynamics of KMP is very similar to that of SSEP: 

V\p,p]e- Na V<d = [ T>[p, p]e- N J ^Vp-W (140) 



where we have introduced the Hamiltonian density: 

n = \ [ P 2 wp 2 - wpwp] (wi) 

To compute the large deviation function T(p* ) , one has to solve the corresponding Hamilton equations 

p=^Ap-W[p 2 Wp}; i5=-±Ap-p(Vp) 2 (142) 

with the boundary conditions 

p(x,0) = p(x) = (l-x)p + xpi) p(x,T)=p*(x) (143) 
p(0,t) =p ; p(l,t)= Pl ; p(0,t)=p(l,t)=0 (144) 

Once again, the last equality simply says that no fluctuations are allowed at the contact with the reservoir. 

8.1. Connection with the SU(1,1) spin chain 

KMP is related to SU(1,1) spin chains in a rather subtle way [19, 20]. Starting from the SU(1,1) coherent 
states for spin k 

\z)= n 1 -, k e zK+ \0) (145) 
(1 — zz) K 

one gets the following expression for the pseudo-spin operators 

7 7 A -4- 7 7 

(z\K+\z) = 2k ; (z\K-\z)=2k ; (z\K z \z) = k— — (146) 

1 — zz 1 — zz 1 — zz 

of the SU(1,1) group. By analogy with equation (40) for the SU(2) case, the Doi-Peliti variables are defined 
through 

<t>* - z; 4> = (147) 
1 — zz 

The difference with the SU(2) case stems from the fact that the energy density variables p, p have to be 
identified directly with the Doi-Peliti variables: 

p = ^=^-_- p = cj>*=z (148) 
1 — zz 

In particular, from (146) and (148), one sees that p corresponds to K~ and not to the z component of the 
spins as was the case for the SSEP and the SU(2) representation [44]. Defining the classical spin 

K x ' y ' x = — (z\K x > v '*\z) (149) 

one then checks that the Hamiltonian (141) corresponds to the continuous pseudo-spin chain 

H = -\ [(Vif x ) 2 + {VK y ) 2 - (VK Z ) 2 } = ~ [VK+ ■ VK~ - (VK Z ) 2 } (150) 
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.2. Remarkable change of variables 



As we did in section 6 for the SSEP, we shall use the symmetry of the evolution operator under simultaneous 
SU(1,1) 'rotation' of all 'spins' to find a basis where a non-local mapping to equilibrium is easily revealed. 
Making a reflexion with respect to the x — z plane maps K y in —K y and lets (150) invariant. As expected, 
it is thus a symmetry of the action. In the (z, z) coordinates it reads 

(151) 

= z! (152) 



z = z: z = z 
By analogy with (147), one defines 
z' 



1-z'z' 1 

The transformation then reads in density variables 

p = <p\i + #*'); p = 



b' = p(l+pp); 4>*' = - 



P 



1 + T ' T 1 + pp (153) 

As for the SSEP, we could now work with these new 'Doi-Peliti' like variables, but to make contact with the 
solution of BDGJL we use slightly different notations 

F = 6*'\ F = 6' 



so that the overall mapping of the action reads 

T 

S[p, p] 



dx 



^-log^ 
F F 



i) 



<b-<l/ <J FF - ^F 2 (VF) 2 + ^VFVf\ 



The fields F, F are related to energy densities through 
P 



F = 



F = p(l + pp) 



so that the spatial boundary conditions (143) are given by 



F(0) = F(l) = 0; 



F(0) = po; 



F(l) = Pi 



(154) 

(155) 
(156) 

(157) 
(158) 



8.3. Non-local mapping back to equilibrium 
Let us introduce the non-local variables 



F' = VF; 
which maps the action into 
S[p, p] -- 



F = V 



F' 



Ax 



F F 



log/" — F'F' 



1 ^.2, 



dxdt^F'F' - ^F r {VF' f 



1 



VF'VF' 



(159) 

(160) 
(161) 



Comparison of (160) and (155) reveals that (159) is a non-local symmetry of the action. The boundary 
conditions (158) now reads 

rL * VF" 

F' = p L -po VF'~-^ = (162) 

As the classical trajectories in the original F, F variables satisfy 
1 



F = |AF + FVF 2 



F = — AF 
2 



F Z VF 



(163) 
(164) 
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we see that on the boundaries, F — and F = implies AF = VF' = 0. Together with the r.h.s. of (162), 
this implies 



VF' = VF' = 



(165) 



and all the currents vanish on the boundaries. (159) thus maps the chain into an isolated one. One can 
continue the mapping p, p — > F,F 



p' = F'(l + F'F'); 
One then gets an overall action 

S[p, p]= dx — 



F' , F' to introduce new energy densities p',p': 
F' 



1 + F'F' 



log - log/" - F'F' -<L+\ g£- 



F 1 



F' 



J dxdt^p'p' -\p' 2 {Vp'f + \vp'V P ' 



(166) 



(167) 



(168) 



which shows that the hydrodynamic limit of the KMP model driven out-of-equilibrium can be mapped back 
to equilibrium through a non-local change of variables. This chain has a detailed balance symmetry, and as 
before its instantons are time-reversal of relaxations and are thus given by 

p' = ~Ap' (169) 



From the equations of motion (142), one sees that this corresponds to 

p 2 Vp = Vp 
In the F',F' variables, this reads 

VF' = 

Mapped back to the initial variables, one gets the instanton equation 

AF 

F = 



(VF) 2 

or, for the density variables, 



P 



F — F 



AF 
VF 2 " 



(170) 
(171) 

(172) 
(173) 



This is the counterpart for KMP of the instanton equation (110) for the SSEP, and corresponds to the 
equation found by BDGJL. The action of this trajectory is 

S[p, p}= J dx [! - log | - log VF] * (174) 

which is precisely the large deviation function obtained in [3] . One sees that once again a non-local mapping 
back to equilibrium enables one to find the instanton equation and thus compute the large deviation 
function. This result can be extended to the whole class of systems defined by a Hamiltonian density 
H[p,p) — a(\7p) 2 /2 — X7p\7p/2 where a is a second order polynomial in p. 
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9. Non-interacting particles in an arbitrary potential driven out-of-equilibrium 



Apart from the simple example addressed in the introduction, the mappings to equilibrium we have presented 
so far only apply at the level of large deviations. We shall show below that such a mapping can also 
be constructed for the full probability distribution, without coarse-graining, of a model of non-interacting 
particles driven out-of-equilibrium. We shall do the mapping in two ways: at the level of probabilities in 
subsection 9.1, and at the level of evolution operators in subsection 9.2. 

We consider an open chain of L sites (index 1 < k < L) in contact with two reservoirs. Each particle at 
site k can jump to a neighboring site k ± 1 with rates W k ^ k ±\. The system is coupled to reservoirs at the 
two boundaries (sites 1 and L) through transition rates Wo->i, Wi_>o, Wl^,l+i, Wl+i^l (See figure 3). 

I 1 1 1 " 1 1 1 1 1 1 1 1 

1 k L 

W^o W L+1 ^ L 

Figure 3. Open chain in contact with reservoirs: generic rates for individual particles. 



A simple case for which the steady state is known is that of equilibrium systems, where probability 
currents vanish. For non-interacting particles, the corresponding detailed balance relation is equivalent to 
balancing the stationary mass fluxes over each bond. Introducing the average occupancies P k q , it reads: 

W ^i = P^Wi^o P c k q W k ^ k+1 = P^ W k+1 ^ k W L ^ L+1 P c L q = W l+ i^l (175) 

The first two equations imply that P^ q takes the form: 

P k q = TT for Kk<L (176) 

e=o "-t- 1 -^ 

Combining the last equation of (175) and (176), we see that this solution is consistent as long as 
W ^ 1 W 1 ^...W L ^ L+1 = 

W L+1 ^ L . . . Wi^o { 1 

This condition can be violated in many physical situations as for instance when a current is forced by 
the reservoirs. In such cases, the steady state occupancies are not the P k q, s and one has to resort to other 
methods for their determination, as has been done for instance by Derrida in [8] for a particle hopping in 
a periodic potential. To our knowledge, the case of open systems with non-vanishing steady current has 
not been considered in the literature and could be interesting in the context of non-equilibrium disordered 
media [5]. 

In the following section, we show how to map our open system back to equilibrium, allowing us to get 
an explicit expression for the steady state distribution. 

9.1. Transformation of probabilities 

Let us consider a probability distribution obtained as a product of Poisson distribution in each site 

L pn k 

P(n 1 ,...,n L ) = l\^ T e- p " (178) 

Its form is preserved by the time evolution and the Pfc's evolve with the conservation equation 

d t P k = -(Jfc+i - J k ) (179) 

where one has introduced the currents 

Jk = -P k W k ^ k -i + Pfe_iW fe _i^ fc (for 2 < k < L) (180) 
J a = -PxWx^o + W ^i J L +i = -W L+1 ^ L + P L W L ^ L+1 (181) 
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In general, the steady state is not given by cancelling all currents J k , as this takes us back to the 
detailed balance conditions (175). In the spirit of previous sections, our aim here is to map our model to 
an isolated system, for which currents vanish at the boundaries. The construction of the equilibrium model 
follows closely the one of the Fokker-Planck equation presented in the introduction. We introduce primed 
occupancies 

h = (Pkiir'Pk+i {prr l p k \ (fon<fc<i-i) 

p> = (prr'Pi - (PoT 1 ; PL = OfT+ir 1 - (pD-'Pl 

where one has extended the definition (176) of P^ q to -P^+i ano - se ^ = 1- Defining the new rates 

WSUk+i = W k+1 ^ k ; W' k+1 ^ k = W k+1 ^ k+2 (183) 
we now have the evolution for a chain of L + 1 sites (0 < k < L) 

d t P' k = -{J'k+i - 4) (184) 
where the currents are defined by 

4 = 4+1=0; and J k = -P k WU k -i+PLiWLi-,k, (forl<fc<L) (185) 
Strikingly, the currents J' and J' L+1 vanish and the primed chain is thus isolated. 

Note that although we mapped our initial open chain of L sites into an isolated chain of L + 1 sites, 
there is no contradiction when counting degrees of freedom. Indeed, from the definition (182), we see that 
the primed occupancies satisfy 

X> fe = (^rM^r 1 ass) 

fc=0 

which means that there are only L independent occupancies in the isolated chain. 

The (equilibrium) steady state of the primed chain is found by imposing J' k — for all fc's: 

P> ocP' (k) TT WU+1 (187) 

Mapping back to the original variables, one obtains the expression of the steady state occupancies P|* 

(188) 



L Z L 



L peq fc-1 peq 



where Zl is a normalization constant given by 

*-£idW (18<J) 

The result (188) is reminiscent of that obtained by Derrida in [8], the differences highlighting the role of the 
reservoirs. 

To highlight the connection with the case treated in introduction, let us consider explicitly particles 
diffusing in a discrete potential, where the rates Wu^>k±i are given by 

In the bulk, they obey detailed balance with respect to the Boltzmann weight P^ q oc e~^ Vk , but the system 
is driven out-of-equilibrium as soon as Vl+i ^ Vq (see relation (177)). The averaged occupancies P k evolve 
with 

pd t P k = Pk^e-i^-^+Pk.xe-^-^-Pk [ e 4m + .-^ +e -|«^-.-^)] (191) 

Taking the continuum limit (k = xL with < x < 1, L 1) and rescaling time with L 2 (diffusive scaling), 
one recovers the Fokker-Planck equation (1) with fi" 1 = T, upon gradient expansion of (191). In particular, 
taking the continuum limit of the microscopic steady state (188), one recovers as expected the result (7) for 
the non-equilibrium Fokker-Planck steady state. Last, with the potential V k defined in (190), the primed 
equilibrium law reads Pc q (fc) oc exp{|(Vfe + Vjt+i)}, which is analogous to the sign change of V for Fokker- 
Planck (see equation (3)), but also induces a smooth averaging of the potential over two neighboring sites. 
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9.2. Transformation for evolution operators 

As shown above, the average occupancies P k fully determine the steady state of non-interacting particles and 
can be computed through a mapping to equilibrium. Such approach does not apply for exclusion processes, 
as exclusion correlates particles. To cast the mapping to equilibrium in a form similar in both cases, we shall 
now work directly with the evolution operator of non-interacting particles. Beyond the quest for a formalism 
ultimately applicable to microscopic models of interacting particles, this point of view proves to be much 
stronger than the one developed in section 9.1, even for non-interacting particles. Indeed, it applies when 
the initial distribution is not factorized, i.e. not in the form (178), and gives access not only to the steady 
state, but also to the whole dynamics. 

For the system of non-interacting particles considered in this section, the probability P{n, t) evolves 
with (we refer to section 2.2 for the notation): 

L-1 

d t P{n,t) = J2{ n k+i W k+i^kP(rik ,n+ +l ) + n+ W k ^ k+1 P{n+ ,n^ +1 ) - [n k W k ^ k+1 
fe=l 



+n k+1 W k+1 ^ k ] P(n) } + W^ [nfP(nf) - mP(n)] + W^i [Pfc) - P(n)] 



+W L ^ L+1 [n+P{n+) - n L P{n)\ + W L+1 ^ L [P(n7j - P{n)] (192) 

Let us introduce the usual Doi-Peliti [13, 32] creation and annihilation operators aj,,a k , defined by 

a\\n k ) = \n k + l); a k \n k ) = n k \n k - 1) (193) 

The master equation can then be written as dt\ip) = — H|t/>), where = J2 n -P( n )l n )> anc ^ the evolution 
operator HI reads 

L-1 

H = y^(<4 + i - a\){a k+ iW k+ i^ k - a k W k ^ k+ i} 
fc=i 

-(4-l)(WW-WW*l) - (a{-l){W L+1 ^ L -W L ^ L+1 a L ) (194) 

Because of the boundary terms, this operator does not correspond to an equilibrium dynamics. We know 
however that its ground state can be mapped to that of an equilibrium operator (see section (9.1)) and it is 
thus quite natural to investigate the question as to whether this mapping extends to the whole operator. 

For the exclusion process, the equilibrium model was constructed at macroscopic coarse-grained level 
through the use of canonical transformations, which mapped the action onto that of an equilibrium, isolated 
system. At the level of operators, canonical transformations correspond to similarity transformations 
H' = Q~ l WQ (see Appendix D). We will now show that using such transformations, one can map HI to 
an equilibrium operator. As these transformations do not modify the spectrum, the determination of the 
eigenstates reduces to the determination of the spectrum of an equilibrium operator. 

We first translate the at's by 1, using the transformation obtained from Qo = e~^fc=i a fc (one has 
Qo 1(1 k Qo = a k + !)• Rearranging the sum, we get for H = Qq 1 M.Q 



L-1 

■£■ 

k=1 



. a{{(W k+ i^ k a k+1 - W k ^ k+ ia k ) - (Wk^k-ia-k - Wk-i^kO-k-i) } 
-a\[{W 2 ^ 1 a 2 -W 1 ^ 2 a 1 ) - (WUoOi - W ^i) } 

-a[{(W L+1 ^ L ~W L ^ L+1 a L ) - (WL^L-iaL-WL-i^ L aL-i)} (195) 

Then, we use the similarity transformation induced by Qi = Y[ k =i (-Pfe q ) 0fc ° fc , which yields 
j Qi X a k Qi = P k q a k 

{Q^a{Qi = (Pk q r^ 



i . (196) 
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and thus for H a = Q^ 1 U Q 1 : 

L-l 

Hi = - al^Wk^k+i (afc+i - a k ) - Wfc-^-i (a k - a fe _i) j 

fc=2 

-4(^2(02-^) - w^ai-ocr 1 )} 

-alj^^+^lP^J- 1 -a L ) - W£_i_n(oi-o £ _i)} (197) 

Noticing that the boundary terms (second and third lines of (197)) look like the bulk term (first line), with 
{PfP') 1 pl a Y m g the role of an operator ao, and (-P^+i) -1 the role of an operator Ul+i, we add two sites 
k = and k = L+l, with their corresponding creation and annihilation operators a^+i- We a l so define 
two vectors \L) and \R), satisfying a \L) = {P^y 1 \L) , a L+1 \R) = (P^) -1 \R) [48]. Defining, on the 
extended space (of L + 2 sites) 

L 

H 2 = -J2 a l{ w k^k+i (ofc+i - Ofc) - VK fc ^ fe _i (a fc - a fe _i)} (198) 
fc=i 

we observe from (197) that 

W 2 (\L)®\ ni ,---,n L )®\R)} = |i)®(H 1 |m > ... > n i ))®|iZ) (199) 

In other words, the action of the extended operator H2 on states of the form \L) ® \rii, . . . , til) (gi \R) reduces 
to that of Hi on the physical state \ni, . . . , til). 

A last transformation brings H2 into an equilibrium form; it is generated by the non-local operator 

ala„ 



Q 2 = cxp y 2™ (200) 

0<p<q<L+l H * 

which transforms the (L + 2) creation and annihilation operators according to 

L + l 

Q 2 X a k Q 2 = y a p ; Q 2 1 a\Q 2 = a\ - a\_ x (l - S k ,o) (201) 

p—k 

These relations are analogous to the non-local primed variables we have introduced for exclusion processes. 
Moreover, the operator H 3 = Q 2 ~ 1 W 2 Q 2 reads 

L-l 

H 3 - y (4+i - a{)[W k+1 ^ k a k+1 - WJLk+i°*} ( 202 ) 

fc=0 

This expression corresponds to an isolated equilibrium process of primed rates W, given in (183). We have 
thus shown that the evolution operator H can be mapped to an equilibrium operator H3 acting on a larger 
space. We shall now explain how to use this mapping to determine the spectrum of H from that of H3. 

The operator H3 acts on a space of L + 2 sites (0 < k < L + 1) but, as seen from its expression, it 
only describes hopping among the first L+l sites (0 < k < L), implying that the last site (k = L + 1) is 
completely isolated. Its eigenstates therefore take the form 

L 

\P*- X ) = (®\P k X ))®\fi + i) (203) 

fe=0 

where f^+i is arbitrary. To determine the eigenstates |P A ) of H, one maps back the eigenvectors of H3 with 
Q 2 and retains only those which satisfies 

Q 2 \P 3 ' X ) = \L) ® |P 2 ' A ) ® \R) (204) 

One then has 

\P X ) =QoQi\P 2 X ) (205) 
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Let us illustrate this procedure to determine the steady-state of H. The degenerate ground state of the 
equilibrium operator H3 takes the form 

|P 3 <°) = ((g)|^)U|/£ +1 > (206) 



k=0 



where \V k ) is a Poisson distribution of mean density [iP k i M being an arbitrary constant: 

, \n k ) (207) 



in.> = E e ^ pr '— q 



Tiki 

rik 



Forcing that |P 3,0 ) satisfies (204) constrains both fi and 

m = L , ; \h+i) = 2^ e 1 L+1> — — - — K+i) (208) 

l^p=Q F k n L + 1 UL+1 - 



One finally checks that applying further Q0Q1 to |P 2: °) as in (205) gives back the steady state obtained in 
section 9.i. 

We have thus shown in this section that the evolution operator of non-interacting particles diffusing on a 
one-dimensional lattice with arbitrary rates can be mapped onto an equilibrium operator. The determination 
of both steady state and excited states can be mapped to an equilibrium problem. The new microscopic 
feature of the transformation (20f ) is that the system has been supplemented with a site at each end, which 
accounts for the additional constant of motions of the isolated system — the total mass is conserved. 

Last, from the knowledge of the mapping to equilibrium for the microscopic dynamics, one can extract 
its counterpart for the hydrodynamic limit. Indeed we show in Appendix D how similarity transformations 
for the evolution operator can be read as canonical changes of variables in the action. From the expressions 
of Qq, Q\ and Q2, one can thus construct a mapping to equilibrium at the level of the action and then takes 
its continuum limit, (see section 2.3.5). 
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10. Conclusions 

The macroscopic fluctuation theory, as developed by BDGJL, consists of two distinct steps. The first is 
the recognition that the coarse-graining level, as measured by the box size, is a parameter playing the role 
of h in a quantum system: the dynamics, including large deviations, becomes in the hydrodynamic limit a 
Hamiltonian dynamics. This step is general, and leaves us with a classical field theory to solve. 

The second step, that has been followed in a number of systems, is to calculate explicitly the trajectories 
starting from the stationary situation and ending in an given configuration at long times (i.e., an excursion 
from equilibrium), and hence obtain the large-deviation function from the action. This can be done trivially in 
systems satisfying detailed balance, just by reversing the corresponding trajectory that relaxes to equilibrium. 
In systems driven out of equilibrium, the obvious detailed balance symmetry is broken. For the driven SSEP 
and the KMP models, however, a hidden detailed balance symmetry can be found explicitly, and used to 
compute the large-deviation functions. In this paper we have shown that the reason for this unexpected 
symmetry, and hence the solvability, is that these systems can be mapped back into their equilibrium 
counterparts. This realization suggested to look back at the simpler case of independent particles diffusing 
in a potential, with sources at the ends. We have found that the same strategy can be applied in this 
elementary case. 

The question that arises is how general this mapping is. The first problem that comes to mind is how 
'one-dimensional' this mechanism is. One can also ask what happens with more general one-dimensional 
cases. A way to start investigating the first question is to consider the simple case of non-interacting 
particles, but in higher dimensions. On the other hand, more general one-dimensional models might perhaps 
admit a mapping back to equilibrium if one is prepared to pay the price of dealing with spatially non-local 
interactions. 
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Appendix A. Coherent state path-integral 

The functional expression for spin operators presents a few mathematical subtleties [17]. Recently Solari [36], 
Kochetov [28], Vieira and Sacramento [39] derived an expression for the action and the associated time 
boundary conditions. We follow here the clear presentation of Stone, Park and Garg [17]. 

Usually, one seeks to calculate the 'propagator' between two normalized coherent states [28] 

(*V r V) (A.l) 
As usual in the construction of a path integral, the time interval [0, T] is divided into N segments and one 
inserts at each time interval L representation of the identity (23), one for each site of the lattice. Letting N 
go to infinity, one ends up with the functional representation of the propagator 

(z f \e- T "\z l ) = [ VzVzaip[-S] (A.2) 

1 Zkik - z k z k 



§= (i + ^ fc (T))(i + ,- fc (o) 4 ) f T dt 

r (i+44)(i+44) J o 



2 1 + ZkZk 



H(z,z) 



(A.3) 



Note that the temporal boundary term differs slightly from those in [28, 17] as the representation we used 
of the SU(2) group is not the usual unitary one (cfr Appendix C). The role of the Hamiltonian is played by 
the quantity 

H(z,z) = -±(z\H\z) (A.4) 

which is computed using (28) or similar expressions for higher powers of the spin operators (see [17]). Zk(t) 
and Zk[t) are two independent complex fields [36]. In the construction of the path integral, one notes that 
•Zfc(O) = zl and ffe(T) = zl [41], while z k (T) and z k (0) are unconstrained [36, 28]. 

Here, we rather wish to compute the physical propagator between two states {n*\ and \ri l ), 
corresponding to fixed initial and final number of particles in each site. The propagator represents the 
probability P(n/, T; n 1 , 0) of observing the system in state (n{, . . . , n£) at time T, starting from (n\ , . . . , n l L ) 
at time 0. Using 2L representations of identity (23), we write 

P{n f ,T;n\Q) = {n f \e- Tk \n l ) = f JJ d/i(z()d/i(4) (n f \z f )(z f \e- TA \z l )(z l \n l ) (A.5) 

J it- 
Keeping in mind that we are aiming to describe a hydrodynamic limit, we introduce the discrete densities 

Pk = n k /2j (A.6) 

One sees by comparing equations (A.2) and (A.5) that the physical propagator is obtained by subtracting 
log(z t \ii l ) + \og(nt \z*) to the action (A.3), and by integrating over z l , z l , , . As we aim to describe the 

large j, L limit, we can use Stirling's formula to obtain the asymptotics of ( 2 \ ), to get: 



P(p',T;p\0) = J J]d/x(zf)d M (4) J VzVzexp{-S} (A.7) 



S = -i5>g[(l + z{z k {T))[l + z k (0)4)] + 2j I dt 



k 



T 



1 z kZk — Z k Zk n ,t- \ 

2^ l + z kZk 



+2jJ2{~Pk log 4 " Pi l °s4 + Pi^gpl + (1 - p{) log(l - p{) + log[(l + z[z[){\ + 44)]} 
k 

The complex fields 4> 4' z k ana - 4 have now to be integrated upon, which can be done by saddle point. 
First differentiating the action (A.7) with respect to z\, z[ yields 

i ds _ 4 pj i ds _ zj P { 

2j dz\ i + 44 4 2j dz { i + 44 4 { ' ' 
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We thus obtain the initial and final conditions 



i + 44 



Pk 



-f i 

44 _ j 



i + 44 



Pk 



(A.9) 



Extremalizing with respect to zj. and z k has to be done carefully as the time integral in the action gives a 
non-zero contribution (See for instance equation 3.9 of [17]). Such a computation leads to 



Zk(0) 



1 OS 



2jdz\ 1 + 5*4 1 + ^(0)4 



1 dS 

2jd4 



z k (T) 



i + 44 i+z k (T)4 



which fixes 



Zfc(0) 



/ 



z k (T) = z[ 



Putting everything together, the action reads 

_ z k z k 
i. ' 



S[ 



z,z;n f ,ni 



n B (z, z) 



i 



L-l 



1 + z k z k 
ZkZk+l {zk+i — Zk) 2 



logz k - log(l + z k z k ) 



ZkZk 



ZkZ k 



(A.10) 
(A.ll) 

H(z, z) 



y — 

2 ^ (1 + z k z k )(\ + z k+1 z k+1 ) \l + z k+1 z k+1 l + z k z k 



Zk+l 



(z k+1 - z fc )(A.12) 



In the context of quantum mechanics, the coherent state (stereographic) parametrization in terms of z k (t) 
and z k (t) is usually transformed into spherical polar coordinates through z k = er % ^ k cot z k = e l ^ k cot 
In our context however, it will prove more convenient to introduce a new parametrization [42] 

Pk 



I- Pk 



Zk 



(A.13) 



Using : f k - Zk = p k we see from (A. 9, A. 11) that the temporal boundary conditions on the fields can be written 

-1 \ %k %k 



as 



p(0) = p i , p(T) = p f , with p(0), p{T) unconstrained (A. 14) 

This highlights the correspondence between the field p and the actual density of the system. The Hamiltonian 
(A. 12) then reads 

H = H B +H +H L , (A.15) 

L — l 

Hb^EK 1 ' W - 1] + Pk(l - Pk+i) [e^-r* - 1] } (A.16) 



k=l 



Ho = a(l - pi)(e pl - 1) + 7 pi(e-^ - 1) 
H L = 6(1 - p L )(e^ - 1) + 0p L (e-P L - 1) 



(A.17) 
(A.18) 



where Hb describes the interaction in the bulk whereas Ho and Hl result from the coupling to the reservoirs. 
Using 



1 z k z k - z k z k . „ r „ 1 -, 

t; t , = PfePfe - PfcM + o log 1 ~ Pk)\ 

2 1 + z fe z fc 2 

we check that all boundary terms cancel so that the classical action reduces to 



S[p,p;Pt,Pf] = 2j 



dt 



^hPk - U{p,p) 



(A.19) 



(A.20) 
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Appendix B. Hydrodynamic limit 



Appendix B.l. Continuous limit of the action 



Let us start from the microscopic Hamiltonian (A. 15) for a system of size L. In the large j limit, each site 
contains a large number of particles, so that the density field tends to self-average. It is thus natural to 
assume the gradients to be small. Expanding the Hamiltonian up to second order in (pk+i — Pk), (Pk+i, Pk), 
one gets 

L-l 



E 

k=l 



^ PfcC 1 — Pk)(pk+i - pk) 2 - 7j(/5fc+i - Pk){pk+i - pk) 



(B.l) 



(B.2) 



One can then introduce a rescaled space variable 
k 

x k = 7T 

2j 

which goes from l/2j to L/2j and becomes continuous in the large j limit. Assuming a diffusive time scale 
and explicitly rescaling the gradients 

L 



1 



Pk+i - Pk^ 77- Vp, p k+1 - p k 
transforms the action into 

Slp(x),p(x)} = (2j) 2 



1 „. 1 



k=l 



L/2j 



dx, dt -> (2j) 2 dt, (B.3) 



dt- 



L/2j 

dx 



pp--a(Vp) 2 



Ho[p, p] = a [1 - p(0)](e m - 1) + 7P(P)(fi~ m ~ 1) 
H x [p, p}=S[l- p(l)](e m - 1) + Pp(l)(e- K1) - 1) 



2jH [p,p}-2jH 1 [p,p}j 
(B.4) 



It is however more convenient to have a continuous variable x going from to 1, so that the extensivity of 
the action is explicit. This can be achieved by a further scaling 

2 



, L 

— x, t — ► — 



2j 



t 



(B.5) 



which maps the action into 

S{p(x),p(x)]=2jLldt 



dx 



PP- 



-cr(Vp) 2 + -VpV/5 



LHofapl-LH^p] } (B 



Some comments are in order. First, the bulk integral is proportional to 2jL, which is consistent with the 
scaling of a large deviation function. Then, the boundary terms scale as 2jL 2 and deserve some further 
analysis. 



Appendix B.2. Spatial boundary conditions 

Let us first note that if the boundary rates a, /?, 7, S are of order 1/L, then the contributions of TIq and TLl to 
the action are of the same order as that of the bulk term TLb- This means that trajectories with fluctuations 
of order 1 at the boundaries give rise to non-vanishing contribution at the level of large deviation, i.e. their 
action scales as 2jL. In particular, such rates would not lead to the spatial boundary conditions (48). 

In the usual case where the rates are of order 1, we shall show below that trajectories with fluctuations 
at the boundaries are forbidden at the level of large deviations. This accounts for the claim in [1, 3] that 
trajectories which do not satisfy strictly the boundary condition correspond to infinite values of the large 
deviation function. 

Let us analyze in details the dynamics at the boundaries. In the large j limit, the probability is 
dominated by trajectories that extremalize the action (see section 4). Using expression (A. 15) for the 
microscopic Hamiltonian, the classical equations read: 
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9(Hi+Hb) 1 - Pi s!-d 2 Pin \ &-pi , ft \ pi -pi m 7 \ 

pi = = — ^ — P2e pi p 2 - — (1 - ^Je^ 2 P1 + a(l - pi)e pi - 7pie P1 (B.7) 

api 2 2 

v d(Hi+Hs) P2< Bi-Sa ii l"~P2 r pjj-p! n , r pi n r -pi n 

= = — le 1 1 — 1 pl — II + a c^ 1 — II — r y[e 11 — 1 

api 2 2 

In the large 2jL limit, one checks that d Pl H.B and d^Ti-B are of order The r.h.s. of equations (B.7) 

are thus dominated by the boundary terms so that, at first order, they read 

pi = a(l - p x )e, h - lPie~ h (B.8) 
' Pl = a^ 1 - 1] - -/[e-P 1 - 1] 

For sake of clarity, we now drop the index '1'. Equations (B.8) can be solved and give 

p(t) = [p (l-p l )(l-e- rt )+p l e-P'(p + (l-p Q )e- rt )] x 

[(1 - Po )(l - e rt ) + e A (p + (1 - p a )e rt )] (B.9) 

m Poe^+l-p + (l-p ) e rt (e^-l) 
e = 1 I& I 1 1 TTfTTT^ il ( BJ0 ) 



^[p] = - lim — logP[p] = o .lim — Straj - oo (B.12) 



where r — a + 7, po = ck/T and pi,/5i are the initial conditions of the fields. For such a trajectory, the 
contribution St of the left boundary to the action is 

S L = -2jtr(l - e"*) [p (l - (H)eh - Pi (l - Po )] (B.ll) 

In the hydrodynamic limit, the time is rescaled by L 2 so that ~ 2jL 2 . We know from section 4.3 that 
the probability to observe a given profile is obtained from the exponential of the action. For non-zero St , it 
is thus of order exp(— 2jL 2 ) and such profile is even more rare than large deviations, whose probability scale 
as exp(— 2jL). Quantitatively, such profile has an infinite large deviation function 

— logPfpl = lim . 

Ijh^roa 2jL Ijh^oc 2jL 

In addition to an infinite large deviation function, the trajectories (B.9) can also present a diverging 
density. Indeed, for large times, one gets 

p(t) ~ e rt po(l ~ Po)(l ~ e-K)[ Pi + (1 - Pl )e p >] (B.13) 

To keep the density finite, one thus needs 

Pi = or pi= 1 (B.14) 

1 _ e p. 

The second solution corresponds to an action St = 2jtT which is once again of order 2jL 2 in the 
hydrodynamic limit. Such a trajectory is thus forbidden at the level of large deviations. We are left with 
Pi = 0, which corresponds to p(t) = po + e~ Tt (pi — po). The corresponding contribution St vanishes exactly 
and in a macroscopic time of order L the density p(t) of the first site gets equal to po = jj^p— • 

The same analysis holds for the right boundary and we have thus shown that for boundary rates of 
order 1, the fields p, p have to satisfy the following boundary conditions at the level of large deviations: 

p(0,t)=p(l,t) = 0; p(0,t)=p = ^—; p(l, t) = Pl = — — (B.15) 

a + 7 + p 

whereas for rates of order 1/L, fluctuations are allowed at the level of large deviations. 
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Appendix C. Unitary and non-unitary representations 



Appendix C.l. Change of basis 

We wish to make the link between the representation of SU(2) matrices we have used (16) and the usual 
"quantum" representation. In terms of the action of these operators, we start from the action of the 
matrices (16) on the occupation kets \n) 

S+\n) = (2j-n)\n + l) (C.l) 
S~|n)=n|n-1) (C.2) 
S z \n) ={n-j)\n) (C.3) 

and perform a similarity transformation to obtain the canonical unitary representation of "quantum" 
operators, that we denote Sf, S*: 



Q-'S+Qln) = S+\n) = ^{n + l)(2j - n)\n + 1) (C.4) 

Q^S-Qln) = S~\n) = y/n(2j - n + l)\n - 1) (C.5) 

Q- 1 S'Q\n) = S;\n) =(n-j)\n) (C.6) 

The usual magnetic number m is related to the occupation number n through m = n — j. We first remark 
that only S + and S~ have to be changed. This suggests we find Q as a function of the number operator 
n = j + S z : 

Q = q(h) (C.7) 
From (C.4-C.5) we have to solve 

= y/n(2j -n+1) (C.8) 



q(n) 



q{n - 1) 



q(n) 

which is done by choosing 



(2j-n)-^-L- = ^/(n + l)(2j-n) (C.9) 
q(n + 1) 



2) ( c - 10 ) 



Appendix C.2. Terms arising from the similarity transformation 

We now denote the basis |n) of kets for occupation numbers as \0), with j cosOk = 2nk — 1 on each site. We 
are interested in the propagator 

P(n',t\n,0) = (n'\e- t "\n) (C.ll) 

where H is expressed in terms of the non-unitary representation (16). However, the field-theoretic 
construction of the action corresponding to (C.ll) is well suited only for 'quantum' spin operators acting on 
kets as in (C.4-C.6). To make the bridge between these representations, let us write 

P(n', t\n, 0) = (n'lQe-^Q-^n) (C.12) 

where H q is now expressed in terms of 'quantum' spin matrices, that is to say, in terms of the matrices 
Sg~' z , which act on \n) according to (C.4-C.6). Inserting two resolution of the identity yields 

P(n',t\n,0) = J di f di f di % diie- s{0 ' fi) (C.13) 

with 
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But Q is diagonal and real- valued on the basis \n) so that 

'2/ 



S) = - fi ' ft + 5 - Hy + S(f " (ai5) 

We thus conclude that using the non-unitary representation of SU (2) yields a new term in the action (the 
ratio of binomials in (C.15)), that we took into account to obtain (A. 3). 

Appendix D. Prom microscopic detailed-balance relation to the symmetry of the action in 
the path integral 

Let us first discuss in detail the relation between similarity transformations of the evolution operator and 
canonical changes of variable in the action, and then show the consequences of the existence of symmetries 
for the system. We then analyze the case of the detailed balance relation. To remain in the context of this 
article, we present below the case of the SSEP in the hydrodynamic limit although the results only rely on 
the existence of the path-integral representation and are thus much more general. 
We consider the propagator between two physical states 

G(p f , Pi ;T) = {p f \e- TH \ Pi ) (D.l) 

where H is the evolution operator. We know from section 2.3 that the path integral representation of this 
propagator is given by 

G(p f , Pi ;T) = Jv[p,p]e-^ LSH ; S H = J dtdx{pp - H[p, p]} (D.2) 

where Tt[p, p] = — -^j(z\H\z) and we used the relations (32) to introduce densities. Let us introduce the 
operator obtained after a similarity transformation 

H = Q^HQ (D.3) 
(D.l) can then be written 

G(p f , Pi] T) = (p f \QQ- 1 e- TH QQ- 1 \ Pl ) = { Pf \Qe- Til Q\p t ) (D.4) 
The path-integral representation of (D.4) leads to 

G( P f, Pi ;T) = J V[p' lP ']e-^ LS »- S„ = _ J^pogfitf + J dtdx{p'p'-H[p'.,p']}{V.$) 

where Ti[p' , p 1 ] = — jj(z\H\z) and Q = (z\Q\z). It is now simple to check that the actions (D.2) and (D.5) 
are related via the canonical changes of variables induced by log Q [21]: 

1 d f 

dx P'p' = TTTTi lo § 2 + / dx PP ( D - 6 ) 



2jL dt 

If Q is a symmetry of the evolution operator, we see that the commutation relation [H, Q] = can 
also be written H = Q~ 1 HQ. At the level of the action, it thus means that log Q induces a canonical 
transformation that leaves the action invariant up to boundary terms. 

The case of detailed balance relation is slightly different as it is not a symmetry of the evolution operator 
but rather a connection between H and its adjoint H\ as we recall. Indeed, for any configuration p 1; p 2 it 
reads 

( Pl | e - tff |p 2 )P cq (p 2 ) = (p 2 \e- tH \ Pl )P cq ( Pl ) (D.7) 
Taking the adjoint of the r.h.s, one gets 

( Pl \e- tH P cq -P cq e- tHi '\p 2 ) (D.8) 
As this holds for all p 1 , p 2 ,t we see by deriving with respect to t and putting t — that 

fft = P-^Peq (D.9) 
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At the level of the actions, it then reads 

J pp-H[ P ,p] = -^L[io g p eq ]{ + j p P -rO[ P ,p] (d.io) 

In the case of the SSEP with periodic boundary conditions, one further knows that 

P cq [p] OC e -VLfd X plo S p+(l-Pn°s(l-p) ( D .H) 

and 

H[ PlP ] = ^aWp 2 -\vpVp- re[p',p'] = \^p> 2 + \v P 'Vp'- (D.12) 

We thus see that the action of H and W are related through the canonical changes of variable induced by 
(D.6): 

p = p; p' = p-—— (D.13) 
I- p 

By further taking a time reversal transformation 

t-vT-t; p'^-p' (D.14) 
one maps back TO to 7i and thus obtain a symmetry of the action, which finally reads 

p m (t)=p(T-t); = -p(T - 1) + 1 P } T ~ * (D.15) 

We get back the transformation (79), as expected. 
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